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f~| We explore the phenomenology of nontrivial quantum effects on low-energy gravity. These ef- 

fects come from the running of the gravitational coupling parameter G and the cosmological 



constant A in the Einstein-Hilbert action, as induced by the Renormalization Group (RG). The 
Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these 
quantum effects, and it is assumed that the dominant dark matter-like effects inside galaxies is 
due to these nontrivial RG effects. Here we present additional details on the RGGR model appli- 
cation, in particular on the Poisson equation extension that defines the effective potential, also we 
re-analyse the ordinary elliptical galaxy NGC 4494 using a slightly different model for its bary- 
onic contribution, and explicit solutions are presented for the running of G and A. The values of 
the NGC 4494 parameters as shown here have a better agreement with the general RGGR picture 
t-H for galaxies, and suggest a larger radial anisotropy than the previously published result. 
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1. Introduction 

In Refs. [1, 2, 3, 4, 5, 6] we presented new results on the application of renormalization group 
(RG) corrections to General Relativity in the astrophysical domain, in particular on a possible 
relation between RG large scale effects and dark matter-like effects in galaxies. The resulting 
phenomenological model was named RGGR. These developments were directly based on the RG 
application to gravity of Ref. [7], and are consistent with the phenomenological consequences 
of diverse approaches to the subject, including the related to the asymptotic safety scenario of 
quantum gravity, see in particular Refs. [8, 9, 10]. 

Currently, in the context of quantum field theory in curved space time, it is impossible to con- 
struct a formal proof that the coupling parameter G is a running parameter in the infrared. However, 
this possibility can not be ruled out. The possibility of General Relativity being modified in the 
far infrared due to the renormalization group (RG) has been considered in different contexts, for 
instance, [11, 12, 13, 14]. The previous attempts to apply this picture to galaxies have considered 
for simplicity point-like galaxies (e.g., [7, 8]). We extended previous considerations by identifying 
a proper renormalization group energy scale pL and by evaluating the consequences considering the 
observational data of disk [1] and elliptical [2] galaxies. We proposed in Ref. [1] the existence of 
a relation between jj. and the local value of the Newtonian potential ( this relation was reinforced 
afterwards [15]). With this choice, the renormalization group-based approach (RGGR) was capable 
to mimic dark matter effects with great precision. Also, it is remarkable that this picture induces a 
very small variation on the gravitational coupling parameter G, namely a variation of about 10~ 7 
of its value across a galaxy (depending on the matter distribution). We call this model RGGR, in 
reference to Renormalization Group corrected General Relativity. 

2. The RGGR effective potential and the running of G and A 

This section is, in part, a review on the RGGR dynamics, but the route to deduce the effec- 
tive potential is new, and an explicit expression for the modified Poisson equation (for spherical 
symmetry) is here presented. Additional details on this deduction will appear in Ref. [16]. Also, 
from this route comes the equation that governs the A running inside galaxies, which was used in 
Ref. [2]. Another approach to find the effective potential, which was used previously, is to use a 
conformal transformation. 

2.1 Effective potential and the modified Poisson equation 

The gravitational coupling parameter G may behave as a true constant in the far IR limit, 
leading to standard General Relativity in such limit. Nevertheless, in the context of QFT in curved 
space time, there is no proof on that. According to Refs. [7, 4], a certain logarithmic running of 
G is a direct consequence of covariance and must hold in all loop orders. Hence the situation is as 
follows: either there is no new gravitational effect induced by the renormalization group in the far 
infrared, or there are such deviations and the gravitational coupling runs as 
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Equation (2.1) leads to the logarithmically varying G(/i) function, 

l + vln(Ai 2 /^ 2 ) 

where /lo is a reference scale introduced such that G(jUo) = Go- The constant Go is the gravitational 
constant as measured in the Solar System. Actually, there is no need to be very precise on where G 
assumes the value of Go, due to the smallness of the variation of G; moreover, within the weak field 
limit, the value of jUo is dynamically irrelevant, as it will be shown. The dimensionless constant v is 
a phenomenological parameter which depends on the details of the quantum theory leading to eq. 
(2.2). Since we have no means to compute the latter from first principles, its value should be fixed 
from observations. The first possibility, namely of no new gravitational effects in the far infrared, 
corresponds to v = 0. 

The form of the action in question is simply the Einstein-Hilbert one, in which G appears 
inside the integral, namely, 1 

c 3 [R-2A , . 

SSA = m J-S-J=i**. (23) 

In the above, G and A are external scalar fields, that is, their dynamics do not come from pure 
classical arguments, but are given from the RG equations; in particular G satisfies (2.2). For the 
problem of internal galaxy kinematics, the cosmological constant effects are negligible (as expected 
and as shown in [2]). 

The above, together with some energy-momentum tensor T^ v , leads to the following field 
equations (see, e.g. [17, 2]) 

, , SkG 

G^i V +Agp V + GDG g^y — GV^VyG =—T ilv . (2.4) 

To find solutions for static spherically symmetric space-times (a natural symmetry for some 
ellipticals, including NGC 4494), the following line element will be considered, 

ds 2 = -e K{r) c 2 dt 2 + e^ r) dr 2 + r 2 dd 2 + r 2 sin(0)J0 2 . (2.5) 

To be clear, by weak field it is mean that both K and tj are much smaller than one. Also it will be 
assumed "small RG corrections", i.e. that G(fl), as given by (2.2), satisfies \G(h)/Gq — 1| <C 1. 

From the geodesies, one finds that K acts as the effective potential, in the sense that its gradient 
is proportional to the acceleration felt by a test particle. But contrary to General Relativity, K* does 
not satisfy a Poisson equation. By solving the field equations in the presence of "dust" of density 
Pif), 

_ SxGo 

.Mo,' ' 



V 2 



, u 
fc + 2vln 1 



c 2 



'-p. (2.6) 



Therefore the relation between the effective potential <J> to the Newtonian one <J>a? can be stated as 

(2.7) 

where jUi is an arbitrary constant (it needs not to be /Xq) and <J> = Kcr/2. 



<& = <S> N -2vc 2 \n ( — 
VMi 



1 We use the ( — h++) space-time signature. 
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2.2 The energy scale setting 

In order to derive a test particle acceleration, we have to specify the proper energy scale p. for 
the problem setting in question, which is a time-independent gravitational phenomena in the weak 
field limit. This is a recent area of exploration of the renormalization group application, where 
the usual procedures for high energy scattering of particles cannot be applied straightforwardly. 
In [1] we introduced a p. identification that seems better justified both from the theoretical and 
observational points of view. The characteristic weak-field gravitational energy scale does not 
comes from the geometric scaling 1/r, but should be found from the Newtonian potential <J>a?, the 
latter is the field that characterises gravity in such limit. Therefore, 

If / would be a complicated function with dependence on diverse constants, that would lead to a 
theory with small (or null) prediction power. The simplest assumption, p oc c^, leads to p oc \Jr 
in the large r limit; which is unsatisfactory on observational grounds (bad Newtonian limit and 
correspondence to the Tully-Fisher law). One way to recover the Newtonian limit is to impose a 
suitable cut-off, but this rough procedure does not solves the Tully-Fisher issues [7]. Another one 
is to use [1] 

PL ( $> N S 



* v*».' ' (29) 

where <I>o and a are constants. The precise value of <I>o is irrelevant for the dynamics in the weak 
field limit, since /Xo the dynamics in this limit does not depend on po as argued above. The relevant 
parameter is a. 

The parameter a is a phenomenological parameter that needs to depend on the mass of the 
system, and it must go to zero when the mass of the system goes to zero. This is necessary to have 
a good Newtonian limit. From the Tully-Fisher law, it is expected to increase monotonically with 
the increase of the mass of disk galaxies. For more details, see [6, 2, 18]. 



2.3 The running of A 

The running of A can be stated as a differential equation that depends on G and R alone. From 
the energy-momentum conservation (T^ v v = 0), together with the Bianchi identities (Guy = 0) [2], 

V v = ^VyG- 1 , (2.10) 

since (DV V — V V D) G 1 = R VK V K G l . The above result is exact. Equation (2.10) also appeared 
in Refs. [19, 20]. And from eqs.(2.9, 2.10), 

<J>'„r 

A' = 2va^ 

Hence, from the knowledge of the matter distribution, one can find <t>w and consequently A(r). 



4nGp 



+ A + 3va 



V 2 <f> N 



(2.11) 
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2.4 Application to spherically symmetric Jeans equation 



Since elliptical galaxies are mainly supported by velocity dispersions (VD), the main equation 
for galaxy kinematics in this case is the following expression for the projected (line-of-sight) VD 
(see also [21]), 

2^0 r ~f r \ t(r)M(r) 



I(R) 



r \RJ 



where 



K(u) = . u 



= V- 1 

2 



3 \ ^rQS-1/2) 



1 1 

2'2 



- dr, 



— B 



>a z 2 2 



(2.12) 



(2.13) 



B(x,a,b) = Jo f 0_1 (l —t b ~ l )dt is the incomplete beta function, T is the Gamma function, r stands 
for the deprojected (spherical) radius, R for the projected (line of sight) radius, l(R) is the luminos- 
ity intensity, £(r) is the luminosity density (found from the deprojection of /(/?)), P is the anisotropy 
parameter (it is zero if the galaxy has an isotropic VD profile) and M(r) is the total (effective) mass 
of the system at the radius r. See Ref. [2] and references therein for additional details. In the case 
of Newtonian gravity without dark matter, M(r) would stand, within a good approximation, as the 
stellar mass inside the radius r, i.e., M* (r). 

For RGGR without dark matter, the total mass inside the radius r is the baryonic mass M*(r). 
Nevertheless, the non-Newtonian gravitational effects of RGGR for spherically symmetric systems 
can be understood from the Newtonian perspective as if the total mass was given by the following 
total effective mass [2] 

M{r)=M Sr (r)+M RGGR (r), (2.14) 



with 



M RGGR (r) 



avc L 



(2.15) 



The above is the the essential equation to model elliptical galaxies (compatible with the spherical 
approximation) within RGGR. 



3. NGC 4494 
3.1 Mass models 

NGC 4494 is an ordinary elliptical recently analysed within RGGR [2], and which has received 
considerable attention previously, in particular due to its apparent lack of dark matter [22]. Here 
we present a variation of the analysis presented in that Ref. [2], namely we do not directly use the 
photometric data with the Sersic extension to model the stellar mass, instead only the Sersic profile 
that best fit the surface brightness of this galaxy is used (an approach here labeled "Full Sersic"). 
See Table 1 and Fig. 1 . 

In part this is relevant to evaluate to what extent our results presented in [2] are robust to small 
changes on the baryonic mass content. The main results are i) within the 1 a error bars, the results 
are compatible; ii) the best fit value for the anisotropy parameter considerably moves towards radial 
anisotropy (from about isotropy to 0.44) (see further comments in the Conclusions). 
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NGC 4494 



RGGR without dark matter 



Stellar model W 


av x 10 7 


P 


T */(||) 


Mass ( !0 10 M e i 


X 2 




Full Sersic+jSp] 


0.77 ±0.6 





4.27±0.31 


10.1 ±1.5 


9.02 


0.38 


Full Sersic+j8p j] 


2 6 +0 " 
z - D -2.o 


44+0-52 


t 4+0.98 




5.10 


0.22 


Full Sersic+K.IMF+j3 [0] 


0.84 ±0.43 





4.23 ±0.21 


10.0± 1.3 


9.09 


0.36 


Full Sersic+K.IMF+j3 [0 A] 


1.17±0.70 


0.19±0.19 


4.09 ±0.32 


9.7± 1.5 


7.06 


0.29 



Table 1: NGC 4494 results for RGGR, assuming that the stellar mass profile is only given by a Sersic profile. This 
table extends a table on this galaxy in Ref.[2], see this reference for further details. (1) j3[ ] indicates isotropic VD, j8[q i] 
indicates constant anisotropy with j8 € [0, 1] (except for the error bars, the above also holds for j8 e [— 1, 1]), K.IMF 
is a reference to Kroupa IMF, and it means that the expected value of Y* was used in the % 2 definition (it was used 
4.2± I.OMq/Lq^ [23]). For details on this procedure see Ref. [2]. 




Figure 1: Two NGC 4494 mass models within RGGR and without dark matter. The curves refer to mass 
models composed by the stellar component (given by a Sersic profile) and RGGR gravity. The black solid 
lines are the resulting VD for each model, the yellow dashed and blue dotted lines are respectively the stellar 
Newtonian and non-Newtonian contributions to the total VD. One of the models assumes isotropy (/3 = 0), 
the second assumes j3 € [—1,1]. The vertical dashed line signs the radius above which the observational data 
is considered for the fitting procedure (10 arcsec). 



The four models presented in Table 1 lead to total stellar masses that are inside the range 
expected from the Kroupa (and Chabrier) IMF's [23], namely M* £ [8.0M Q , 12.0M©]. The found 
anisotropy values are also reasonable [23]. On the other hand, only the second model has a value 
of a compatible with the baryonic matter amount of this galaxy 2 , all the others display a "dearth of 
RG effects" for the given baryonic content. This has some similarities with the ACDM case, where 
some radial anisotropy is needed to find reasonable amounts of dark matter [22, 24]. 

3.2 The running of G and A inside NGC 4494 

From the NGC 4494 results above, here the corresponding running of G and A inside this 
2 In accordance with expectations driven by the analysis of other galaxies, see the conclusions. 
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galaxy are shown. A priori, according to eqs. (2.2, 2.9), G(r) depends on the value of <I>o (or 
Ho), nevertheless dynamically such constant plays no role (it can be anything and it will not have 
dynamical impact whenever the weak field approximation holds). 




Figure 2: The plots show the variation of G (left) and A (right) inside NGC 4494, both plots considers 
two models, "Full Sersic+j3p]" (black on the left plot, blue on the right plot) and "Full Sersic+j3[Q,i]" (cyan 
on the left plot, black on the right plot). Both of the G plots in the left are normalised with G(0) = Go; 
the latter fixes the (dynamically irrelevant) constant <t>o for each model. The right plot compares the ratio 
of the absolute value of A(r) and 4-nGp^/c 2 to A(25/? e ) (for further details, see also [2]). The solid curve 
corresponds to 4-nGp t (r)/(c 2 A(,), the dashed curve to A(r)/Ao. In both cases, the absolute value of A, at 
any radius, is too small to significantly change the average internal galaxy dynamics, which is in accordance 
with the approximations used to derive the effective potential of RGGR. 



4. Conclusions 

In summary, RGGR is a model based on the theoretical possibility that the beta function of 
the gravitational coupling parameter G may not be zero in the far infrared. Currently, there is no 
way to directly deduce this behaviour from first principles, nevertheless eq. (2.1) is a natural, if not 
unique, possibility that has appeared many times before in the context of Quantum Field Theory in 
curved space-time. 

RGGR without dark matter is a model with one phenomenological free parameter (a) which 
is capable of dealing with the kinematics of diverse galaxies. The a relation to other physical 
parameters is being disclosed in a work in progress [18]. In particular, we are unveiling a correlation 
between a and the galaxy baryonic mass. The galaxy that was in worst agreement with such 
correlation was NGC 4494, since it displayed a too low av value for its mass (considering the 
solutions presented in Ref.[2]). Here it is commented that with a small change on the modelling of 
the baryonic mass (i.e., modelling its mass exclusively from its Sersic profile, which is well inside 
the eiTors on the conversion from light to mass), a greater tendency towards higher radial anisotropy, 
lower mass and higher a was found (model with free constant anisotropy). The corresponding mass 
is still inside the mass range associated with the Kroupa IMF [23], and fits much better within the 
general picture that is being unveiled [18]. 
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